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ABSTRACT 

Context. Long timeseries of data increase the frequency resolution in the power spectrum. This allows for resolving stochastically 
excited modes with long mode lifetimes, as well as features that are close together in frequency. The CoRoT fields observed during 
the initial run and second long run partly overlap, and stars in this overlapping field observed in both runs are used to create timeseries 
with a longer timespan than available from the individual runs. 

Aims. We aim to measure the mode lifetimes of red giants and compare them with theoretical predictions. We also investigate the 
dependence of the mode lifetimes on frequency and the degree of the oscillation modes. 

Methods. We perform simulations to investigate the influence of the gap in the data between the initial and second long run, the total 
length of the run and the signal-to-noise ratio on the measured mode lifetime. This provides us with a correction factor to apply to the 
mode lifetimes measured from a maximum likelihood fit to the oscillation frequencies. 

Results. We find that the length of the timeseries, the signal-to-noise ratio and possible gaps do have a non-negligible effect on the 
measurements of the mode lifetime of stochastically excited oscillation modes, but that we can correct for it. For the four stars for 
which we can perform a fit of the oscillation frequencies, we find that the mode lifetimes depend on frequency and on degree of the 
mode, which is in quantitative agreement with theoretical predictions. 

Key words, stars: late-type - stars: oscillations - methods: observational - techniques: photometric 
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QQ 1. Introduction the evolutionary phase of the star. This frequency and mode de- 

fvj pendence of the mode lifetime has not yet been observed in red 

LlJ The uninterupted long timeseries of high-precision photomet- gj^^f^ ^jjj investigated here. 

. ric data obtamed with the CoRoT satellite (Baglin et al. 2006!) The mode lifetimes (r, e-folding time for the amplitude) can 

opened many new possibilities for red-giant asteroseismol- measured from the width (FWHM) of the frequency peaks 

^ ogy. Red giants have extended turbulent atmospheres in which ^„ oscillation mode in the power spectrum as 
^—^ solar-like oscillations are stochastically excited and intrinsically 

I \ damped. These oscillations have periods of a few hours, and 1 

^ . can have mode lifet imes ranging from days to hundreds of days ^ ~ "nf" 

' dPupret et al.ll2009h . The mode lifetime of a stochastically ex- 
?S ' cited mode is the characteristic time it takes for an excited mode To measure the mode lifetime it is therefore important to have 
^ to damp exponentially. This is an important parameter as it is sufficient frequency resolution to resolve the linewidth of the os- 
■ ■ direcdy related to the mode inertia and therefore the cumu- ciUation mode. Examples of stars with long (> 50 day) and short 
lated work integral. From this work integral we can derive in (< 20 day) lifetimes were first shown by De Ridder et al. (2009D 
which part of the star the oscillations a re driven and d ampe d (see and Carrier et al. (2010) and investigated for a larger number of 
iDupretet al.ll2009l for more details). iDupret et all (12009 ) also stars by ,Baudin et al., (201^). Where .Baudin et alj dZOK)) pro- 
predict that mode lifetimes depend on the degree of the mode Pose an automated procedure to fit the line width and height for 
and on frequency. This dependence follows similar trends for *e three most prominent oscillation modes in a star, with the 
stars in different evolutionary phases, although the exact varia- Purpose of a statistical study of the amplitudes of the modes, we 
tion of the mode lifetimes with frequency changes depending on adopt here an individual analysis for a smaller sample of stars, 

but with longer timeseries. The data we have at our disposal are 

* The CoRoT space mission which was developed and is oper- described in Sect. 2. Because we do not have 'infinitely' long 

ated by the French space agency CNES, with participation of ESAs timeseries and there are considerable gaps in those series, we 

RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, perform simulations to identify a possible bias on F in Eq. [T] 

and Spain. Light curves can be retrieved from the CoRoT archive: This bias will be taken into account in computing the 'unbiased' 



|http://idoc-coro t.ias.u-psud.fr, mode lifetime. These simulations are described in Sect. 3. 
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We performed a detailed analysis of four stars. For each star 
we then computed an echelle diagram (frequency versus fre- 
quency modulo large frequency separation, Grec et al. 1983]). 
From this echelle diagram we aim to identify the degree of the 
oscillation modes, or at least identify modes with the same de- 
gree to investigate the dependence of the mode lifetime on de- 
gree. This is possible when the observed oscillation modes have 
high- order (») and lo w-degree (0 and follow the asymptotic re- 
lation 1^11^1 ([l980|), for which we use here the following equa- 
tion: 



v„,t^ Av\n + y + e\-£(£+l)Do, 



(2) 



with V the oscillation frequency, Av the large separation between 
modes of same degree and consecutive orders, which is inversely 
proportional to the sound travel time through the star. Do is sen- 
sitive to deeper layers in the star and e to the surface layers. The 
large separation can be estimated from the relation between the 
frequency of maximum oscillation power (v^ax) and Av as found 
in an ensem ble study of the first long run of CoRoT observations 
(lHekkeretalJl2009l) : 



0.78 



Av - Av, 



(3) 



or sim ilar relations by IStello et aP (l2009t) and iMosser et alJ 
( 1201 0). We use this relation to check whether we indeed found 
Av and not a (sub)multiple of Av. 

In general three small separations are defined: 6vo2 is the 
spacing between adjacent modes with { = Q and 1-2, Sv\-i is 
the spacing between adjacent modes with ( - \ and ^ = 3, and 
6vq\ is the amount by which the £ - \ modes are offset from the 
midpoints between the { - Q mode on either side. If Eq.|2]holds 
then Sv{)2 - 6Do, 6v\3 = IQDq and 5voi - 2Do, although there 
is ev idence that fo r the majority of red-giant stars 6voi is nega- 
tive (ICarrier etal.f 2010) and Bedding et al. (2010), i.e. the ^ = 1 
modes are located on the right side of the midpoints between the 
{ - modes instead of on the left side of the midpoint as is the 
case for the Sun. This is most likely caused by stellar evolution 
and the asymptotic relation given above may not be valid for 
evolved stars. 




Fig. 1. Colour-magnitude diagram of the field stars observed dur- 
ing LRaOl using the 2MASS / and Ks photometric passbands 
(black dots) (Deleuil et al. 2009). The blue diamonds represent 
the red-giant stars also observed during IR, the green triangles 
represent the red-giant candidates which show clear power ex- 
cess and the red crosses represent the four stars for which we 
could fit the individual oscillation modes. 
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2. Observations 

For the present analyses we use CoRoT data taken during the 
Initial Run (IR) and first long run in the anti-centre direction 
(LRaOl) in the exo-planet field. For information on CoRoT and 
technical details about the CoRoT data reduction we refer to 
iBagUn & Chaintreuill (1200^ . The IR lasted 58 days from 2"^^ 
February to March 2007 during which the satellite was 
pointed towards a field with the centre at equatorial coordinates 
(a, S) = (102.60°, -1.7°). During the 131 daylong LRaOl, which 
lasted from October 2007 to March 2008, the satellite observed 
a field around equatorial coordinates (a, 6) = (102.72°, -0.2°). 
The typical timestep of the observations is 32 s, although some 
targets are observed at 512 s intervals. 

Visual and near-IR photometry are avai lable for all stars in 
the IR and LRaOl field (iDeleuil et al.ll2()09h. These colo i irs are 
affected by reddening, but as shown by 'Bess ell & BrettI (Il988l) 
near-IR colours are least affected and provide a first estimate of 
the spectral type. A J - K versus K colour magnitude diagram is 
shown in Fig.[T] 



Fig. 2. Window functions for a 394 day timeseries (black solid 
line), the same timeseries with a large gap, i.e., no data between 
day 58 and day 263 (red dashed line) and for the same timeseries 
with the data of the first 58 and last 131 days 'merged' (green 
dot-dashed line). The window functions are oversampled by a 
factor of ten for visual purposes. 



We use these colours to distinguish between main-sequence 
and red-giant stars. In this way we found 305 red-giant candi- 
dates based on colours (blue diamonds in Fig.[T]), which are ob- 
served during both the IR and LRaOl. We computed a combined 
timeseries for these candidates, in which we 'merged' the time- 
series. Simulations to investigate the effect of merging the IR 
and LRaOl datasets are described in Sect. 3. Then we computed 
power spectra and global fits to them. These fits are composed of 
a power law to account for the background of the spectra due to 
effects such as granulation, and a Gaussian fit to t he oscillation 
power excess, similar to the fitting described by iHekker et al.l 
(l2009h . For 24 stars ( green triangles in Fig. [TJ we found clear 
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Fig. 3. Average of the lifetimes observed for 1000 simulated re- 
alisations as a function of their input lifetime. The 'full', 'gap' 
and 'merge' cases are indicated with red asterisks, blue trian- 
gles and green diamonds respectively. Errors on the mean are 
plotted, but are in most cases smaller than the symbol size. The 
black dotted line is the one-to-one relation. 
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Fig. 4. Same as Fig. [3] but here the red asterisks, blue triangles 
and green diamonds indicate the 'full', 'long' and 'short' cases, 
respectively. 



power excess and we investigated the oscillations in more detail. 
For four stars the noise level was low enough to perform fitting 
of individual oscillation modes. These four stars are indicated 
with red crosses in the colour-magnitude diagram in Fig. [T] and 
the detailed analysis of these stars is described in Sect. 4. 

3. Simulations 

The large gap of 205 days between the initial run (58 days) and 
the first long run in the anti-centre direction (131 days) has a 
large influence on the window function (Fig. |2]i by introducing 
prominent sidelobes, and these have to be taken into account 
in the analysis of the data. To investigate how we should treat 
the gap and what effect the gap has on the resulting parameters 
obtained from a fit to the power spectru m we performed sim - 
ulations using the simulator described by IChaplin eta n dl997h . 



First we performed simulations of noise-free timeseries which 
cover the full timespan of 394 days with 512 second cadence. 
Each timeseries has one stochastically excited oscillation mode 
at a frequency in the frequency range of red giants observed with 
CoRoT, i.e., 5, 15, 25, 95 //Hz and for each frequency we 
simulated ten different mode lifetimes of 5 up to 95 days. For 
each combination of frequency and mode lifetime we simulated 
1000 realisations. 

In these simulations of the full timespan we introduced a gap 
by removing all data between day 58 and day 263, but not alter- 
ing any of the time stamps. In addition, we also merge the time- 
series data of the first 58 days and last 131 days by removing 
the gap and altering the timestamps (see Fig. |2] for the window 
functions of these scenarios). The latter is an option for red gi- 
ants as these have stochastically excited modes with finite mode 
lifetimes. These stochastic modes do not have a coherent phase 
and one may therefore in principle merge data sets with different 
epochs. 

We analysed for all three cases, i.e., 'full', 'gap' and 
'merge', the individual realisations. In all instances we fitted the 
power spectrum usin g a maximum likelihood technique (e.g., 
[Anderson et alj|1990l) . where the final model M used for com- 
parison is a convolution of the model with the power spectrum 
of the observed window function of the data set, normalised to 
unit total area (see Eq.|4]l. This takes the redistribution of power 
caused by gaps in the data into account. 



M 



((Vce„,-V)/B,)2j 



* window 



(4) 



In this model, each of the j oscillation peaks are fitted with the 
amplitude (A), the central frequency (Vcen) and FWHM (B) as 
free parameters. To maximise the possibility of obtaining the 
best values, i.e., the solution of the global instead of local max- 
imum, we use different input parameters for the linewidth and 
height of the frequency peak. From the total of twelve fits per- 
formed for each timeseries we choose the one with the maxi- 
mum likelihood. Before we used the result we also tested if we 
could consider the fitted function as a good approximation of the 
parent distribution by means of reduced figure-of-merit function 
(MERITred, see I Anderson et al.ll 19901 and references therein): 



MERIT, 



O i - Mi 

Mi 



(5) 



with Oi the observed power at frequency / and M, the model at 
the same frequency and / the number of degrees of freedom, i.e., 
the difference between the number of independent data points in 
the fit and the number of free parameters in the model. A fit with 
MERITied less than one can be considered a good fit to the data. 

The results immediately showed that the frequency at which 
the oscillation mode occurred did not have any influence on the 
fitting and the resulting lifetime. Therefore we computed the life- 
time of each realisation at one frequency (35 //Hz) and computed 
the mean of these results. The standard error on the mean of the 
sample is computed as crj y/N with cr the standard deviation of 
the sample and the number of measured lifetimes over which 
the mean is computed. These are shown in Fig. [3] This figure 
clearly shows that the measured lifetimes are equal to or lower 
than the lifetimes that were used for the simulations and that the 
offset from the one-to-one relations is larger for the 'gap' and 
'merge' cases than for the 'full' case. For the 'full' timeseries 
the results become less reliable above ~ 60 days while the 'gap' 
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Fig, 5. Same as Fig.[3]for the 'gap' (triangles) and 'merge' (di- 
amonds) cases with 1.5, 1, 0.5, 0.25 and 0.125 times the gap in 
the observations. These are indicated in red, black, green, blue 
and yellow, respectively. The results for the full timeseries are 
indicated with the light blue asterisks. 
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Fig. 6. Top: same as Fig. [3] for the 'merge' case with height-to- 
background ratio (HI B) of the order of 1 0, 25 and 50 indicated in 
red, green and blue, respectively. The black solid line indicates 
the polynomial used to obtain the 'unbiased' lifetime. Bottom: 
the difference between the 'observed' lifetimes and the fit. 



and 'merge' results are consistent only for lifetimes below 20 
days. 

The limited timespan of the observations compared to the 
mode lifetimes could be the cause of the discrepancy between 
our results and the input values. To investigate this further, we 
repeated the simulations and analysis for timeseries with half 
('short') and twice ('long') the timespan of the initial simula- 
tions, i.e., 197 and 788 days and without a gap. We compared the 
resulting average lifetimes computed from the analysis of indi- 



vidual realisations. These are shown in Fig.|4] Indeed we see that 
for the 'long' timeseries the results are consistent with or overes- 
timating the input for all mode lifetimes investigated here, while 
for the 'short' timeseries the discrepancy between input and ob- 
served lifetime start to emerge at lifetimes of about 25 days. 
This confirms that the total time span of the observations {T) 
influences the measured mode lifetime (r), and that we can only 
measure 'unbiased' mode lifetimes when (roughly) T/t > 10. 
We note that the 'gap', the 'merge' and the 'short' timeseries, 
which have similar effective timespans, all show similar results. 

The results for the 'gap' and 'merge' timeseries are approxi- 
mately equal for lifetimes up to ~ 35 days. For higher input life- 
times the mean of the observed lifetimes for the 'merge' time- 
series are systematically lower than for the 'gap' timeseries. The 
reason for this increasing difference is probably due to the differ- 
ent frequency resolutions and might therefore vary as a function 
of the (fractional) length of the gap. 

To test the influence of the length of the gap on the resulting 
lifetimes, we repeated the analyses as described above for time- 
series with a total length of 394 days, and different gaps of 1.5, 
0.5, 0.25 and 0.125 times the length of the gap in the observed 
data. The results are shown in Fig.|5] From these simulations it 
becomes clear that the deviations of the results shown for the 
'gap' and 'merge' timeseries increase with increasing length of 
the gap, i.e. they depend on the effective timespan of the data. 
We see that for timeseries with a shorter effective timespan we 
can measure only shorter 'unbiased' mode lifetimes, which is 
in agreement with the results from the 'short' and 'long' time- 
series. Furthermore, the length of the gap in the observed time- 
series (205 days) is longer than the lifetimes considered here. 
This changes when we investigate the same lifetimes for data 
with shorter gaps. Nevertheless, we find smooth trends between 
the observed and input lifetimes for all gaps and mode lifetimes 
considered here. 

We also investigated the influence of the height-to- 
background ratio of the oscillation modes on the measurement of 
the mode lifetime. We did this with simulations of 'merge' time- 
series, similar to the observed timeseries we analysed, on which 
we added noise such that the height-to-background ratio (H/B) 
of the oscillation mode is of the order of 10, 25 and 50. In these 
simulations the noise has been added to the simulated timeseries, 
while the height of the oscillation mode (H) was derived from 
the root-mean-square of the flux (A) and the input mode lifetime 
(t) as H - A'T. Results for these simulations with different H/B 
are shown in Fig.|6] These results showed that the results are to 
a large extent insensitive to H/B, although the inclusion of noise 
in the simulations provided different results compared to results 
from noise-free simulations (see Fig. [3]l. These differences are 
most likely due to the noise and due to the addition of an extra 
free parameter (background level) in the fitting. 

From these simulations we conclude that we need to apply 
a correction to the 'observed' mode lifetime to obtain the 'unbi- 
ased' lifetime whether we leave the gap in the data or we merge 
the timeseries of the initial and long run. To avoid selection of 
frequencies which have high power due to the more severe side- 
lobes in the window function of the data with a gap, we have 
performed the analysis on observed power spectra of merged 
timeseries. The lifetimes of the observed modes then need to be 
coiTected. For this correction we use a third order polynomial fit 
to the lifetimes measured from the simulations including noise: 



y^5A + Q.96x - Q.QUx^ + 1.1- 10 



(6) 



with y the observed lifetime and x the input lifetime as shown 
in Fig. |6] This shows that the fit is not exact, but that the devi- 



4 



S. Hekker et al.: Oscillation mode lifetimes of red giants observed during the initial and first anticentre long run of CoRoT 



ation from the 'observed' lifetimes is less than five days for all 
lifetimes and HjB. In the analysis of the observed stars we will 
use this correction as a look-up table to compute 'unbiased' life- 
times for observed lifetimes between 5.4 and ~ 35 days, i.e. the 
range of observed lifetimes over which the correction is com- 
puted. For any observed lifetimes outside this range we can infer 
only lower or upper limits. 

4. Individual stars 

From the red-giant candidates selected from near-IR colours and 
power excess as described in Sect. 2 we have four candidates for 
which we can fit the individual oscillations. The fitting is per- 
formed to the non-oversampled power spectra in a similar way 
as for the simulations (Sect. 3). The fitting is performed in two 
steps. First the background, for which we use the power law fit to 
account for effects such as activity and granulation as described 
in Sect. 2 is fitted. Then, we compute a global maximum likeli- 
hood fit to all oscillation modes. In this fit the previously deter- 
mined background is kept fixed, while the frequency, height and 
line width of each oscillation mode and the noise level are free 
parameters. No dependence of parameters on the degree of the 
mode and thus on the azimuthal order is included in the fitting. 
This is both because we do not know the degree of the modes 
and because we assume that the stars are slow rotators and the 
rotational splitting is negligible. 

The selection of the oscillation modes is based on a statisti- 
cal test of the binned power spectrum. For this test we bin the 
power spectrum over intervals of three frequency bins. Then we 
compute the probability of the power in the binned power spec- 
trum to be due to noise. This probability is computed using ?lx^ 
distribution in which we take the width and number of the bins 
into account in the degrees of freedom. Frequencies at which the 
probability of the power not being due to noise is larger than 
80% are selected as candidate oscillation frequencies. After per- 
forming the fit, we also verify that in the ratio of the observed to 
the fitted spectra no prominent peaks are left. Therefore we com- 
pute the relative height {s) for the investigated frequency range 
for which the probability of observing at least one spike with 
this height due to noise is 10%, following the formulation by 
IChaplin et al.l (|2002|) and references therein. 

From the resulting height (//) and width (F) for each oscil- 
lation frequency we also compute the root-mean-squared am- 
plitude (A = VttHT). This parameter is a measure of the total 
power of a mode and reflects the balance between the damping 
and the excitation of the mode. H and F are not independent 
and this is taken into account in the computation of the error in 
A using a co rrelation coeff i cient of -0.9 following the method 
described bv lChaplin et al ] (l2000h . 

4.1. CoRoT 102732890 

CoRoT 102732890 is a star with frequencies showing maximum 
oscillation power at about 26 fiHz. Our statistical test indicated 
nine frequency intervals to have more than 80% probability to 
be due to signal. We fit for these frequencies and the results are 
shown in Fig. |7] together with the ratio of the observed power 
spectrum to the fit. Indeed, in the fitted region no clear oscilla- 
tion modes are present anymore. The MERITied is also below 
one (0.79), which indicates a good fit. However, one could ar- 
gue that the power around 22 yuHz is not due to one oscillation 
mode but due to two modes. Therefore, we also perform a fit 
with ten frequencies. According to the MERJTred this is also a 
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Fig. 8. EcheUe diagram of CoRoT 102732890, using ten fre- 
quencies in black symbols, with a size proportional to the height 
of the fitted oscillations. When the two modes at ~ 22 //Hz are 
replaced by one mode then the two lower left black symbols are 
replaced by the red symbol. 



good fit (0.79). The maximum likelihood is (slightly) higher for 
the fit with ten frequencies compared to the fit with nine frequen- 
cies. We perform a likelihood ratio test (e.g., Appourchaux et sl\ 
1995) to see if the increase in the maximum likelihood justifies 
the addition of three free parameters. In this particular case we 
see that there is a 35% probability to get this increase without the 
addition of an oscillation mode. From this probability we cannot 
draw any firm conclusions and therefore we continue the analy- 
sis using both nine and ten frequencies. The values used for the 
fits are listed in Table [1] together with the computed amplitude 
and 'unbiased' mode lifetime. 

From the observed Vmax and Eq. [3] we expect a large separa- 
tion Av of about 3.2 juHz. Indeed when we plot the ten frequen- 
cies in an echelle diagram with a folding frequency of 3. 15 //Hz, 
the dilTerent mode sequences became visible (see Fig.[8]l. 

Inspection of Fig.[8]shows that the modes are located in two 
roughly vertical regions. The interpretation of this is not unam- 
biguous. One could argue that the broad ridge consists of two 
ridges, namely the { - 2 and £ = ridge, which seems to be the 
case for the fit with ten frequencies (black symbols in Fig.[8]|. We 
could also interpret the broad ridge as the £ = 1 ridge where we 
see more scatter due to le ss eff icient trapping of th e modes, (see 
e.g.. lBeddingetalJ (l2010l) . andlDupret et al.l (i2009h for more de- 
tails on mode trapping). The latter might be the case for the fit 
with nine frequencies as the two bottom black symbols of the 
broad ridge are then replaced by the red point. 

4.2. CoRoT 102788308 

CoRoT 102788308 has a strikingly clear and regular pattern 
around a frequency of about 50 /iHz. As seen in Fig. l9land[T0l 
for the power spectrum and echelle diagram, respectively, and 
Table |2] for the fitted parameters. The MERJTred of the global 
fit to the power spectrum is 1.2, which indicates that this is a 
reasonable fit to the data. When looking at the ratio of the obser- 
vations to the fit we indeed still see some structure, with a main 
feature at about 46.3 //Hz. At this frequency we see for many 
CoRoT stars the fourth harmonic of a one day feature. Therefore 
we do not fit this structure in this case, although we cannot ex- 
clude that part of this signal is due to stellar oscillations. In addi- 
tion, the peaks at frequencies about 35, 40, 65 and 68 //Hz might 
be due to oscillations, but the height to background ratio is not 
large enough to confirm that and therefore these are not taken 
into account in the fit. 



5 



S. Hekker et al.: Oscillation mode lifetimes of red giants observed during the initial and first anticentre long run of CoRoT 

Table 1. Parameters for the fit to CoRoT 102732890. The result between '()' is the oscillation mode at ~ 22 //Hz in the fit with 
nine frequencies, while the results for the interpretation of two oscillation modes at ~ 22 juHz are indicated in square brackets. The 
frequencies, linewidth and height are fitted parameters, the height-to-background ratio (H/B) is measured as the observed height of 
the oscillation mode in the power spectrum over the fitted background. The amplitudes and 'unbiased' lifetimes are derived from 
the fitted parameters as described in the text. 



frequency 


linewidth 


hei^ 


;ht 


H/B 


amplitude 


lifetime 






10^ ppm^//Hz ' 




ppm 


days 


20.96 ± 0.033 


0.026 ± 0.003 


0.5 ± 


0.1 


9.7 


64± 26 


>100 


(22.2 ± 0.2 


0.27 ± 0.04 


0.20 ± 


0.02 


11.6 


130 ±33 


10.0!3«) 


[22.0 ±0.1 


0.040 ± 0.004 


0.37 ± 


0.05 


11.4 


69 ± 17 


>100] 


[22.3 ±0.1 


0.16 ±0.03 


0.22 ± 


0.04 


11.7 


103 ± 30 


28.3+f/] 


24.08 ± 0.06 


0.07 ± 0.01 


0.6 ± 


0.1 


19.5 


107 ± 38 


>100 


25.14 ±0.03 


0.036 ± 0.003 


0.7 ± 


0.1 


30.9 


87 ± 27 


>100 


25.62 ± 0.04 


0.050 ± 0.002 


0.93 ± 


0.05 


30.8 


120 ± 20 


>100 


27.27 ± 0.07 


0.21 ±0.01 


0.3 ± 


0.1 


20.3 


136 ± 83 


16.1!!;^ 


28.7 ± 0.6 


0.2 ±0.1 


0.11 ± 


0.05 


9.9 


85 ± 37 


17 7 +82.3 
^'■'-9.6 


30.45 ± 0.07 


0.079 ± 0.005 


0.14 ± 


0.01 


8.4 


59 ± 11 


>100 


31.6 ± 0.1 


0.084 ± 0.008 


0.22 ± 


0.03 


14.1 


76 ± 18 


>100 




Fig. 7. Top: power spectrum for CoRoT 102732890 (black solid line), with a fit for nine oscillation modes (blue dash-dot line) and a 
fit with ten oscillation modes (red dashed line). Bottom: ratio between the observed power spectrum and the fit with nine frequencies 
and ten frequencies (same colours and linestyles as in the top panel). The value of s for this frequency interval is 8.1. See text in 
Section 4 for a description of s. 



From the Vmax of about 50 //Hz and Eq. |3] we expect a Av 
of about 5.4 //Hz. Indeed we find Av - 5.3 //Hz by folding the 
frequencies. The echelle diagram shows two very clear ridges, 
i.e., the { - and f = 1 ridges (although the actual degree of 
the ridges is not identified) with one point outside the central 
ridge in Fig. (TO] This mode has most likely been shifted due to 
an avoided crossing and in that case the centre ridge is the i = I 
ridge. 

4.3. CoRoT 102762620 

CoRoT 102762620 shows oscillation frequencies at slightly 
lower frequencies than CoRoT 102732890. With the statistical 
test of the smoothed power spectrum we find oscillations roughly 
between 17 and 27 //Hz (see Fig. [TTIand Table[3]for the fitting 
parameters). 

From the frequency of maximum oscillation power of about 
22 //Hz we predict a Av of about 2.8 //Hz. From the observations 



we find a consistent value of 2.7 //Hz for this star. This results in 
one narrow ridge and a second ridge which consist either of two 
ridges relatively close together or one ridge with larger scatter. 



4.4. CoRoT 102767771 

CoRoT 102767771 shows oscillations at about the same fre- 
quencies as 102732890, although the power of the oscillations in 
this star is less than for 102732890. We find a similar value for 
Av of 3.1 //Hz (Fig.IH) as for 102732890, which is in fine with 
the predictions. The echelle diagram shows two clear ridges and 
only the mode with lowest frequency seems to be located outside 
the right ridge. From the cuiTent analysis we cannot say whether 
this ridge is curved and this mode is located on the ridge, whether 
the frequency is altered due to avoided crossing, or whether this 
is induced by noise. 
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Fig. 9. Top: power spectrum for CoRoT 102788308 (black solid line), with a fit for nine oscillation modes (red dashed line). Bottom: 
ratio between the observed power spectrum and the fit. The dashed horizontal line indicates s (see text for more details). 



Table 2. Parameters for the fit to CoRoT 102788308. The frequencies, linewidth and height are fitted parameters, the height-to- 
background ratio (H/ B) is measured as the observed height of the oscillation mode in the power spectrum over the fitted background. 
The amplitudes and 'unbiased' lifetimes are derived from the fitted parameters as described in the text. 



frequency 


linewidth 


height 


H/B 


amplitude 


lifetime 


yuHz 




10^ ppm^Hz"' 




ppm 


days 


42.3 ±0.1 


0.15 ± 0.02 


0.048 ± 0.007 


11.0 


46 ± 13 


33.2!-' 


45.00 ± 0.06 


0.083 ± 0.004 


0.132 ±0.009 


24.2 


59 ± 10 


>100 


47.43 ± 0.06 


0.175 ± 0.009 


0.096 ± 0.006 


32.0 


73 ± 13 


23.2+2-^ 


50.21 ± 0.04 


0.097 ± 0.005 


0.18 ±0.01 


55.6 


74 ± 13 


>100^ 


52.66 ± 0.06 


0.137 ± 0.007 


0.106 ±0.007 


32.7 


67 ± 12 


44 4+13.9 


55.5 ±0.1 


0.17 ± 0.01 


0.040 ± 0.004 


20.9 


47 ± 10 


24-7+3.6 

i+9?.l 


57.2 ± 0.2 


0.3 ± 0.2 


0.023 ± 0.002 


9.1 


49 ±36 


60.8 ±0.1 


0.27 ± 0.02 


0.017 ± 0.001 


9.5 


38 ±8 


10.0+1^ 


63.33 ± 0.08 


0.21 ± 0.04 


0.01 ±0.01 


10.9 


33 ±23 


16.1!!-^ 
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Fig. 11. Top: power spectrum for CoRoT 102762620 (black solid line), with a fit for nine oscillation modes (red dashed line). 
Bottom: ratio between the observed power spectrum and the fit. The value of s for this frequency interval is 8.0 (horizontal dashed 
line). 
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Table 3. Same as Table |2]but for parameters of the fit to CoRoT 102762620. 



frequency 


linewidth 


height 


H/B 


amplitude 


lifetime 


yuHz 




10^ ppm^yuHz"' 




ppm 


days 


17.7 ± 0.2 


0.14 ±0.03 


0.7 ±0.1 


4.5 


54 ± 17 


40 9+^' ' 


18.9 + 0.2 


0.13 ± 0.01 


0.24 ± 0.03 


17.6 


98 ± 24 




20.37 ± 0.07 


0.041 ± 0.004 


1.4 ±0.1 


50.6 


136 ± 30 


>100 


21.7 ± 0.2 


0.13 ±0.01 


0.37 ± 0.03 


47.2 


124 ± 26 




23.1 ±0.1 


0.3 ± 0.3 


0.10 ±0.06 


19.8 


89 ±77 


58.3!38-, 


24.0 ±0.1 


0.13 ±0.06 


0.10 ±0.03 


12.5 


63 ±27 


24.40 ± 0.07 


0.07 ± 0.01 


0.13 ±0.02 


12.8 


53 ± 13 


>100 


25.8 ± 0.6 


0.241 ± 0.005 


0.047 ± 0.006 


8.0 


59 ±21 


12 4+"-^ 


26.9 ± 0.3 


0.38 ± 0.06 


0.05 ± 0.03 


14.8 


79 ±57 


4 8+'-^ 




25 

frequency [/iHz] 



Fig. 13. Top: power spectrum for CoRoT 102767771 (black solid line), with a fit for eight oscillation modes (red dashed line). 
Bottom: ratio between the observed power spectrum and the fit. The value of s for this frequency interval is 8.1. 

Table 4. Same as Table|2]but for parameters of the fit to CoRoT 102767771. 



frequency 


linewidth 


height 


H/B 


amplitude 


lifetime 






10' ppm-yuHz"' 




ppm 


days 


17.23 ± 0.09 


0.088 ± 0.02 


0.12 ±0.05 


5.1 


59 ±25 


>100 


21.31 ±0.09 


0.22 ± 0.02 


0.10 ±0.01 


12.7 


82 ±24 


14.8i» 


22.86 ± 0.09 


0.04 ± 0.01 


0.2 ±0.1 


9.1 


50 ±24 


>100 


24.21 ± 0.09 


1.1 ±0.2 


0.03 ± 0.01 


12.6 


108 ± 40 


<5 


26.0 ±0.1 


0.11 ±0.01 


0.13 ±0.02 


13.2 


66 ± 19 


98 6^'°° 


27.2 ± 0.3 


0.44 ± 0.06 


0.08 ± 0.01 


21.5 


105 ± 26 


Q Q + I.6 


29.1 ±0.2 


0.18 ±0.03 


0.10 ±0.02 


10.9 


74 ±23 


21 9+1'-'* 
2 0+''-'^ 


30.3 ± 0.3 


0.51 ±0.07 


0.034 ± 0.008 


6.1 


74 ±26 



5. Discussion and conclusions 

5.1. Simulations 

Simulations of long (a few hundred days) timeseries of data, 
with different timespans and different length of gaps reveal that 
for short mode lifetimes we can recover the input lifetime of 
the mode, irrespective of the presence of a gap in the data, or 
whether we merge the data. For longer lifetimes we see a general 
under estimation of the 'observed' mode lifetimes compared to 
the input lifetimes, in which the 'observed' mode lifetime seems 
to reach a maximum value. The investigations of longer and 
shorter timeseries reveal that for roughly t < O.IT the observed 
lifetimes are consistent with the 'unbiased' lifetime. For longer 



lifetimes a correction has to be applied, although the uncertainty 
in the resulting values increases with increasing lifetimes. We 
also investigated the trends in the observed lifetimes as a func- 
tion of the length of the gap and whether merging the data has 
a large influence. As expected the shorter the gap the better 
Furthermore, the results for the data with a gap and merged to- 
gether differ for longer lifetimes, but as a correction is necessary 
in both cases, we concluded that we could merge the data of the 
initial run and second long run of CoRoT observations. This has 
the advantage that sidelobes of the window function are avoided 



8 



S. Hekker et al.: Oscillation mode lifetimes of red giants observed during the initial and first anticentre long run of CoRoT 
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Fig. 10. Echelle diagram of CoRoT 102788308, with symbol 
sizes proportional to the height of the fitted oscillations. The 
dashed lines indicate the same ridge with the dashed-dotted line 
the midpoint between them. 




1.5 2.0 
frequency mod 2.7 [yuHz] 



3.0 



Fig. 12. Echelle diagram of CoRoT 102762620, with symbol 
sizes proportional to the height of the fitted oscillations. 
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Fig. 14. Echelle diagram of CoRoT 102767771, with symbol 
sizes proportional to the height of the fitted oscillations. 
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Fig. 15. Linewidth of the oscillation modes of the four stars stud- 
ied in detail as a function of frequency. The dashed lines indi- 
cated linear fits to the results. 
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Fig. 16. Mode lifetimes for ridge 'A' and 'B' in the echelle dia- 
grams of the four stars separately. The degree of the ridges are 
not known and the assignment of 'A' and 'B' to a ridge is ar- 
bitrary. The symbols are the same as in Fig. [15] The oscillation 
modes for which we could only find a lower limit are plotted 
with a lifetime around 100 days, but with a small offset. An off- 
set is also introduced in two of the three modes with a lifetime 
of58.3 days (102762620). 



5.2. Observations 

For four red-giant stars we have data from the initial and sec- 
ond long run of CoRoT at our disposal with high enough signal- 
to-noise ratio to perform a global fit to the individual oscilla- 
tion frequencies. The oscillation modes were identified using a 
statistical test of the binned power spectrum. In this global fit 
the frequency, linewidth and height are free parameters for each 
identified oscillation mode. From the parameters obtained from 
the fitting we also compute the root-mean-squared amplitude and 
the mode lifetime. The latter includes a correction using the third 
order polynomial obtained from the averaged mode lifetimes of 
the simulations including noise of merged data. 



5.2.1. { = 2 modes 

Interestingly for none of the four stars did we find a clear { - 2 
ridge in the echelle diagram. For 102732890 and 102762620 one 
ridge is broad and could be interpreted as an ^ = and £ = 2 ridge 
but for 102788308 and 102767771 only two ridges, separated by 
~0.5Av, are present. The non-detection of the ^ = 2 modes could 
possibly be due to a lack of frequency resolution, in which case 
the { = and f - 2 modes are interpreted as one oscillation or 
due to the lack of height in the power spectrum. We investigate 
both possibilities. 

From the individual oscillations observed for 102788308, we 
can obtain the differences between consecutive frequencies with 
reasonable accuracy. We find here two typical differences: 2.4 
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Table 5. Global parameters for the four stars: frequency of maxi- 
mum oscillation power (vmax) and the large frequency separation 
(Av in yuHz, the ratio of the mean height of ridge 'A' to ridge 'B' 
for which the degrees are unknown, the effective temperature 
(Teff) in K and the asteroseismic mass and radius in M© and Rq, 
respectively. EiTors are typically of the order of 1 and 0.05 yuHz 
for Vmax and Av, 100 K for T^s and below 20 and 10 % for M 
and R, respectively. 



star 


^max 


Av 






M 


R 




//Hz 


yUHZ 




K 


Me 


Ro 


102732890 


25.5 


3.15 


1.1± 0.1 


4140 


1.2 


13.0 


102788308 


51.8 


5.3 


1.6 ±0.1 


4260 


1.3 


9.4 


102762620 


21.9 


2.7 


0.3 ± 0.03 


4355 


1.5 


15.6 


102767771 


25.8 


3.1 


0.4 + 0.1 


4250 


1.4 


13.7 



and 2.8 yuHz. This means that the amount of offset of the £ -I 
modes from the midpoints between the { - Q modes ((5voi) is of 
the order of 0.2 //Hz. If we interpret the centre ridge as f = 1 
ridge (Fig. [TOb then this ridge is mostly located on the left side 
of the midpoint (dot-dashed vertical line in Fig.fTOb of the £ - Q 
modes (dashed vertical lines in Fig.fTOll. The position of the £ -\ 
ridge on the left side of the midpoint between ( = modes might 
indicate that this star is in an earlier evolutionary state and that 
it seems to be justified to use the prediction of the small sepa- 
ration from the asymptotic approximation (Eq.|2|. When we use 
this approximation and the estimate of ^vqi of the observed fre- 
quencies, we find that (5vo2 should be of the order 0.6 //Hz. The 
frequency at ~57.2 //Hz lies by this amount outside the centre 
ridge but no other convincing features are present, although the 
frequency resolution of the dataset is sufficient to resolve this 
separation. 

Also for 102767771 we find evidence for only two ridges 
and we seem to miss the £ - 2 ridge. In this case this might 
be caused by a combination of a lower signal-to-noise ratio and 
the low frequencies of the oscillations. At these low frequencies 
we expect also a smaller separation between the { - Q and £ - 2 
modes compared to the 0.6 //Hz estimated for 102788308 and the 
frequency resolution of the spectrum is therefore more critical 
to detect them, although the frequency resolution of ~0.06 //Hz 
should be sufficient. If the modes are very wide, i.e., short mode 
lifetimes or unresolved mode pairs, a single observed power ex- 
cess hump might actually consist of two modes. Fitting the two 
modes as one will overestimate the line width considerably and 
the correction to compute the 'unbiased' mode lifetime (Eq.|6]l 
might also no longer be valid, as this is determined for a single 
oscillation mode. 

The theoretical computations bv lDupret et alJ (l2009t) predict 
that the trapping of ^ = 2 modes is more efficient than the trap- 
ping of the £ - \ modes due to the increased size of the evanes- 
cent region for modes of higher degree. Due to this more ef- 
ficient trapping the £ - 2 modes are predicted to have similar 
heights in the power spectrum as radial modes. In this predic- 
tion the visibility effects are not taken into account. For the Sun 
it is known that modes with degrees £ - 0, 1, 2 have relative 
heights 1.0, 1.5 and 0.5. These relative heights depend on limb 
darkening and might be different for red giants compared to the 
Sun. For the stars investigated here we computed the mean of 
the fitted heights per ridge and the ratio of these mean values. 
The results are listed in Table |5] and cover a considerable range 
(note that the inverse of the ratio should also be considered as the 
degree of the ridges 'A' and 'B' is unknown). This might indi- 
cate that the relative visibilities of modes with different degrees 



in red giants are different from that of the Sun and this might 
cause the non-detection of E - 2 modes, although this would be 
in contrad iction with re d giants for wh ich £ -2 mod es ha ve bee n 
observed lDe Ridder et alJ (e.g.. l2009l) : ICarrier et alJ (e.g.. l2010l) . 

It might also be possible that the lifetimes of the £ -2 modes 
are long and that these modes are not resolved. Unresolved non- 
radial modes with long lifetimes have smaller heights than the 
closest radial mode (Du pret et al . 2009). With longer observ- 
ing runs, these modes can be resolved, which will improve the 
chances to detect them. 

Despite the difficulties to detect £ - 2 modes encountered for 
the stars investigated here, these modes have been observed in a 
considerable number of (but not all) red-giant stars usin g both 
CoRoT and Kepler data with shorter tirnespan s (De Ridd er et alJ 
'2009':'Carrier et al.'2O10 HBedding et al.l2O10HHuber et al.l2010l: 
Kallinger et al. 2010a). The reason why the £ = 2 modes are not 
observable in the stars investigated here is not entirely clear The 
fact that for some stars £ - 2 modes are observed and that for 
other stars, such as the ones investigated here we cannot observe 
these modes, might also contain information about their internal 
structures. The small number statistics of only four stars does 
not allow us to draw any firm conclusions at this stage. 

5.2.2. Mode lifetimes 

Theory predicts that mode lifetimes depend on frequency and 
on the degree of the mode. In general the mode lifetimes de- 
crease for modes with increasing frequencies and are longer for 
modes with higher degrees (Dupret et al. 2009). We test this pre- 
diction using the four stars we analysed in detail, by plotting the 
linewidth of each mode as a function of frequency and perform a 
linear fit (see Fig.fTSl). For three of the four stars there is a clear 
increase of the linewidth with frequency present, i.e., the mode 
lifetimes decrease with increasing frequency. 

Also a considerable number of lifetimes observed in the four 
stars are also in agreement with the (r) ~ T'^^ scaling predic- 
tion presented by Chaplin et al. (2009). For the effective temper- 
atures mentioned in Table |5] lifetimes of the order 11-15 days 
can be expected. We indeed see in Fig.[T6]that many results are 
present in this regime. 

Furthermore, we investigate the dependence of lifetimes on 
the degree of the modes (see Fig.[T6b. The most obvious exam- 
ple of this dependence is present in 102767771 where all modes 
in the left ('A') ridge in the echelle diagram in Fig.[T4]have life- 
times shorter than 20 days, while all modes in the right ('B') 
ridge have lifetimes longer than 20 days. Although we could 
not identify the degree of the ridges 'A' and 'B', this result is 
consistent with predictions by Dupret et al. (2009), namely that 
modes with different degrees can have systematically different 
mode lifetimes. For the other stars the dependence of the mode 
lifetime on degree is not so clear, although the two longest mode 
lifetimes for 102788308 are both present in the ridge indicated 
with the dashed vertical line in Fig. [10] (ridge 'A'). The three 
modes with the same lifetime of 58 days (102762620) are also 
located in the same ridge. 

From the echelle diagrams we were not yet able to identify 
the degrees of the modes. Because Dupret et al. (2009) predict 
that non-radial modes have longer lifetimes than radial modes, 
we use that as additional information to identify the degree of the 
modes. For 102767771 this would indicate that the right ridge in 
Fig. E] (ridge 'B') is the ^' = 1 ridge and the left ridge ('A') the 
£ - Q ridge. 

The investigation of mode lifetimes as performed here on 
four CoRoT stars will benefit from longer uninterrupted time- 
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series on a larger number of stars. This will be possible in the 
near future. The Kepler satellite will be taking data of a pres- 
elected set of stars during its entire lifetime of nominally 3.5 
years. This increased timespan with only short interruptions and 
the lower noise level of the photometric timeseries compared to 
CoRoT or any other facility cuiTently in operation, will provide 
data which will be even better suited for investigating the mode 
lifetimes of solar-like oscillations in red-giant stars. Analysis of 
the currently available Kepler data with a timespan of ~230 days 
is in progress. 

5.2.3. Stellar parameters 

For the four stars we investigated in detail we also compute 
the asteroseismic mas s and r adius using Av, v,nax and T^ff, as 
described by Kallinger et al.l (l2010b). For this determination 
we need the effective temperatures. The effective temperatures 
have been derived using Ex o-Dat phot ometric data (iDeleuil et al.l 
2009) u sing the rela tions of lAlonsoet al. ( 1999), as described by 
Baudin et al.l ("2010). Because of large, systematic discrepancies 
between the reddening values along t he various lines of sigh t es- 
timated from the ext i nction maps of iDobashi et"an (l2005h and 
iRowles & Froebriciil (l2009h a representative value. Ay - 0.6 
mag, was adopted for all stars. In view of the quite uncertain 
reddening, and although we note that there is a good agreement 
between the values derived from optical and near-IR colour in- 
dices, only the temperatures based on 2MASS data are consid- 
ered for the determination of the mass and radius. The values of 
these stellar parameters are listed in Table |5] 

5.2.4. Comparison with models iDupret et al.l <2009ft 

The frequency of maximum oscillation power provides an 
indication of the evolutionary state of the stars. The stars 
102732890 and 102767771 have a v,™x of about 26 juHz. Stars 
with Vmax roughly between 25 and 40 fjHz are commonly ob- 
served iHekker et al. 2009) and are to a l arge ext ent low-mass 
He-burning stars in the red clump (Mig lio et al.l [2009). Both 
102732890 and 102767771 are therefore'most likely to be red- 
clump stars. The oscillations in 102762620 are ce ntred around 
22 yuH z, which is comparable with model E of Dupret et alj 
(120091). This is model of a 2 star high in the red giant 
branch. For 102788308 we find oscillations with Vmax ~ 50 juHz. 
These fr equencies are comparable to models B, D and E of 
dPupret et al. 2009), which represent a 2 Mq star ascending the 
giant branch, and 3 M© H-shell and He-shell burning models re- 
spectively. The models all have slightly higher masses than the 
asteroseismic masses derived from the observations. The exact 
positions of the stars in the H-R diagram is as yet unknown, 
due to the difference in mass between the models and the ob- 
servations, although we can still conclude that these four stars 
are likely to be in different evolutionary phases on the red-giant 
branch. 

The frequency pattern of 102788308 is strikingly regular and 
the linewidths and trend of the linewidth with frequency for this 
star are different from that of the others. This is most likely due 
to the different evolutionary state of this star This will be further 
investigated using detailed modelling of these stars, which will 
be presented in subsequent publications. 
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